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I. INTRODUCTION 


An analysis of the performance of a local/mobile radio communications 
system utilizing the slotted ALOHA random access protocol is investigated in this 


thesis. 


A. ALOHA NETWORKS 


1. Pure Aloha 

In Aloha networks, the information comes in packets. In a conventional 
pure ALOHA network the users transmit their packets whenever new data are ready 
for transmission. When two or more users transmit their packets simultaneously, a 
collision will occur in the channel, and the packets may overlap partially or 
completely (Figure 1). The collision destroys the information in the packets, and the 
affected users have to retransmit their packets. These collisions cause a major 
degradation in performance, which limits the maximum channel throughput of pure 


ALOHA to 1/2e or 0.184 [Ref. 1]. 


2. Slotted Aloha 
In a slotted ALOHA network, time is segmented into slots. Each slot has 
a duration equal to the packet transmission time or length. Users are required to 
synchronize the transmission of their packets so that the leading edges of the 


packets are aligned with the beginning of the predetermined time slot at the receiver 





Figure 1. The Collision in Pure ALOHA Network. 


Since packets arrival at the receiver are synchronized, when more than one packet 
Is transmitted simultaneously, a collision will occur in the channel, and the packets 
will overlap completely (Figure 2). The elimination of partial packet overlap 
substantially improves the performance of the network. The maximum channel 
throughput is increased to 1/e or 0.368 (twice the maximum channel throughput of 


pure Aloha) [Ref. 1]. This maximum channel throughput is not necessarily the 


tagged packet 





Figure 2. The Collision in Slotted ALOHA Network. 
throughput attained by the system. In heavy traffic, the likelihood of collisions 
increases Significantly. After each collision, the affected users have to retransmit their 
collided packets. As the number of retransmitted packets becomes large, the channel 
will reach a point where a collision occurs every time slot. Hence, no packets can 


be transmitted successfully and the system is said to be unstable. 


B. ALOHA NETWORK IN LOCAL/MOBILE COMMUNICATIONS 

The Aloha network is receiving a great deal of interest for use in applications 
which involve large physical areas (eg. ships, multi-story buildings, warehouses, urban 
environments, etc.) and require random multiple access. The two primary limitations 
of the Aloha network are its poor maximum throughput and its inherent instability. 


However, the Aloha protocol is still desirable because it is easy to implement. 


C. PACKET CAPTURE EFFECT 


1. Advantages 

In local/mobile radio communications, each transmitter produces a 
different signal power level at the receiver due to the differences in receiver- 
transmitter range, shadow fading, and Rayleigh fading. The arriving packet with the 
highest power level will have the best ee to capture the receiver. As a result, it 
is possible for the receiver to accurately decode one of the packets involved in a 
collision. Since the collision does not necessarily destroy all involved packets, the 
channel throughput will increase considerably. The capture effect significantly 


improves the performance of the Aloha network [Ref. 2]. 


2. Limitations 
The advantage of improving channel throughput is achieved in favor of 
users near the base station and at the expense of users far from the base station. The 


user far from the base station will have a low probability of capturing the base 


station receiver, and consequently will have to transmit its packet more often than 


a uSer near the base station and hence suffer a long packet delay [Ref. 2]. 


D. PERFORMANCE MEASUREMENT 

Several methods have been considered to measure the performance of the 
Aloha system. They are packet reception probability, the correct reception of the 
packet address header, and capture measures based on a threshold. Since the packet 
reception probability provides the best measure, this method is selected to measure 
the performance of the system. A transmitted packet is considered to have captured 
the base station receiver only when the entire packet is received correctly in the 


presence of one or more interfering packets [Ref. 3]. 


E. ORGANIZATION OF THE THESIS 

The analysis is divided into four sections. Section II investigates the capture 
effect for a ring model in local/mobile communications. Section III presents the 
numerical results for a three-ring network. Section IV analyzes the system throughput 
for a multiple ring Aloha network using a Markov chain model. Finally, discussions 


and conclusions are presented in Section V. 


II. PACKET CAPTURE PROBABILITY OF A RING MODEL 


The analysis of the packet capture effect in a ring model for local /mobile 
radio systems will be presented in this section. The probability of correct packet 
capture will be evaluated for a system using ideal coherent binary phase-shift-keying 


(BPSK) modulation. 


A. SYSTEM MODEL 

A network consisting of uniformly distributed users in an annular ring will 
be studied. A uniformly distributed density function for the user traffic is 
assumed. The model illustrated in Figure 3 is the system model for this study. 


The traffic generated within a distance r of the base station is 


r 


G(r) = 2n [ xe(x) dx (1) 


0 


where r is the distance from the user to the base stati , The user traffic density is 
g(r). 
The total offered traffic in the network is obtained as G = G(r > %). 


Therefore 


G = 2n if xg(x)dx (2) 
0 


user 


Base station 





Figure 3. System Model 


The spatial distribution function for packet generation ( the probability of a 


a packet generated within a distance R < r) is [Ref. 3] 


F(r) = Priksr} = = xg(x)dx (3) 
0 


Therefore the probability density function for packet generation is found to be 


dF .(r) am 2 


dr 'G 





tio - rg(r) (4) 


The probability density function for packet generation in ring R,, defined as 
T;<rSr, 1, 1S 
FA) 
rSTST,, 


FRATIR) =) F,(r,,,)-F Ar) | 
0 , otherwise. 


(5) 


Assuming that users are uniformly distributed in an annular ring of inner radius 
R, and outer radius R,, we have the user traffic density 
G 
———,Rs<rs<R, 
g(r)=} n(Rj-R) (6) 
0, otherwise 
Substituting (6) into (4), we have the probability density function for packet 
generation as follows 
ar >? R,srsR, 


FAT) Re a 


0 , otherwise 


(7) 


The distribution of packet generation of a ring with R, < r < R, is shown in 
Figure 4. The probability of a packet being generated within a ring of radius R < 1; 


is given by 





Figure 4. User Traffic Density. 


2 G 
F(r) = 2 {x2 a (8) 
‘ G4 (R2-R?) 


which can be evaluated to yield 


2 
r; 








F (7) 7 er (9) 
R,-R, 
The probability of a packet generated within a ring of radius R < r,,, is 
Tie 
F7,,,) = — (10) 


Substituting (10) into (5) and (7), we get the probability density function for 


packet generation in ring R. 


= pissy (11) 


Tia 77; 





f RT [R;) = 


where 1 equals 0, 1, 2,...,.n-1 and R. = {r,<r<vr.,;} 


B. BPSK BIT ERROR PROBABILITY 

We assume that users communicate with the base station using coherent binary 
phase shift keying (BPSK). For ideal coherent binary phase shift keying (BPSK) 
modulation in a Rayleigh fading channel, the received signal of a particular user 0 
in the presence of n-1 interfering users is given by (assuming the background 
Gaussian noise is much smaller than the interference) 


n-l 
r(t) = ¢,8,t) + Vie, e s(t) (12) 


i=] 
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where s(t) are BPSK signals, the c,’s are Rayleigh distributed signal amplitudes, and 
the 8.’s are the uniformly distributed phases of the interfering signals (we set 6) = 
Q and assume that 0,, i=1, 2,..., n-1 are relative to this reference). 


At the sampling time the decision variable of BPSK receiver is 


n-1 
Vp = C4) + >, a; cos8, (13) 
i=l 


where a, = {+1} and are independent and identical distributed random variables 
(1.1.d.). The second term on the right side of the above expression is the sum of 
independent Gaussian random variables. 


The received signal power is 


S = E ((c, a,)*} = E {cg} (14) 
The received noise power is 
N = E {(vy- cy a) (15) 


Substituting equation (13) into equation (15), we have the received noise power 


as follows 


n-1 
N = E ((}\c, a, cos0,)’} (16) 


i=] 


Noting that the mean of cos’6, is 1/2, we find 
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Nie : E {c?) (17) 
i=] 
Hence, the received signal to noise ratio is 
2E {c?} 
SNR = = eS (18) 


The bit error probability of user 0 in the presence of n-1 interfering users for 


BPSK in a Rayleigh fading channel is [Ref. 4] 


1 - ———— (19) 
2 
SNR 


1 
P. = — 
eee) 


1+ 





The near/far effect is modeled by assuming that the received power w conforms 


to an inverse a power propagation law; that is 


w=cr* (20) 
where r is the distance of the user from the base station and a is the propagation 
constant. Generally, 3<a<5 for ultra high frequen ,UHF) applications. The 
constant c depends on the gains, heights, and power of the transmitter and receiver 
antennas. Without loss of generality, we assume that c is normalized to unity for all 
users; that is, users are taken to be identical and use omnidirectional antennas with 
radiation pattern maxima in the horizontal plane. 


From (20), the received power of user i is 


Mee 


w, = E{e;} = 7," (21) 


Substituting (21) into (18), we have 





2w 2r.” 
Fj dpe ne ee (22) 
n-1 n-1 n-] r.\7% 
-¢ t 
Ww, Tr — 
i=1 i=l i=1 * : 
and the bit error probability is 
Pa (23) 
an?) 
which can be written as 
P, = > igo ___ (24) 
n-1 r.\% 
1+y >} 
xe 


For numerical analysis purposes, we select the propagation constant a = 4. 


Hence, 


1 - —————_—_ (25) 


The bit error probability for the case of M-level diversity is [Ref. 3] 


re 


m-1 
M-1+kVl+p\* 
Pi, = PHY ( a ee (26) 


where the value of pu is given by 


> i (27) 





C. CODING FOR PACKET CAPTURE 

Since all packets are transmitted with equal power and because of the near/far 
effect and the Rayleigh fading channel, the operating signal-to-noise ratio can be very 
low. To overcome this problem, we propose a combination of diversity and forward 
error correction. The Viterbi decoding scheme is employed. This scheme proposes 
maximum likelihood decoding, which is relatively easy to implement for codes with 
small memory orders. A binary convolutional code with rate R = 1/2 and constraint 
length 7, and Viterbi hard decision decoding is considered. Define the input 
modulation as M-ary and the output modulation as Q-ary. When binary coding is 
used, the modulator has only binary input (M = 2). Similarly, when binary 
demodulator output quantization is used (Q = 2), the decoder has only binary 
output. In this case, the demodulator is said to make hard decisions (M = Q). To 
simplify the implementation, most coded digital communication systems use binary 
coding with hard-decision decoding. Since errors at the output of the Viterbi decoder 


occur in burst even if the errors at the input of the decoder are independent, the 
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probability of correct packet reception cannot be evaluated by assuming independent 
errors. It has been shown that probability of correct packet reception Pc is lower 


bounded by [Ref. 3] 


Pees (iP) (28) 
where L is the packet length and Py is the union bound first event error probability 


[Ref. 5]. For hard decision decoding, a tight upper bound on P, is given by 
isc: S vapey (29) 
d=d, 


where a, is the number of code words of Hamming weight d, d, is the free distance 
of the convolutional code, and P, is the first event error probability of a code word 


of weight d. The first event error is made with probability [Ref. 5] 


“ (d 
f | Ph Ge as , d odd 
e=(d+ € 
.= , 1)/2 ; 3 (30) 
dj2 _ df2 é = d-e 
— EL 1-P + P 1-P , d even 
lee) 7 Pa «(2 PP 


The bound of (29) can be further simplified by noting that for d odd 


se A 


P,< f 
d 
e=(d+1)f2 \ © 


Bee (iP)? (31) 


Expressions involving P,,, can be moved out of the summation since Pyy is 


independent of the summation index e. As a result, the upper bound for P, is 


ls) 


d 
d 
P, Phy (-P, ee f (32) 
e=(d+1)y2 \ € 


Taking the lower limit of e at 0 instead of (d+1)/2, we obtain 


d 
PreeP (apne > (3 (33) 
e=0 


The sum can now be evaluated to obtain an upper bound on P, of 


P, < 24 Po (1-P,,)” (34) 
Similarly, for d even an upper bound on P, is found to be 
lf d dj2 d\ .an 
fas EO) Pip -P,)” Ge» 
2 a/2 e=d/2+1 € 
This expression can be simplified in a manner similar to that used to obtain 


(34) to get an upper for P, when d is even 


Pied Pega P a2 (36) 
which is identical to (34). 
The correct packet capture probability q(n) is obtain by averaging the 
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probability of correct packet reception over all possibic iocation of n users. Thus 


RR, RF 
q(n) = iF if - [ PRC) <n eee ange (37) 
RR, R, 


where f,(r) is given by (7). Now substituting (30) into (42), we obtain a lower bound 


for the correct packet capture probability as 
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RR, 
1G) 2 | Coa AR AAR ee ee (38) 


R, R, 
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Ii. NUMERICAL RESULTS 


A numerical analysis has been performed. Programs using Matlab provide 
numerical results for the above correct packet capture probability for user 0. Eight 
interfering packets have been used in the analysis. As the number of the interfering 
packets increases the computation time becomes very lengthy. It took approximately 
24 hours to compute q(1), q(2)...q(6). Therefore, the simplified P, in (36) is used to 
calculate q(1), q(2),... q(7), q(8). The result for q(n) using the simplified bound of P, 
in (36) yields results that are very close to those obtained with (30). For each 
additional user after eight interferers, the results do not change significantly. The 
source codes in Appendix A are executed to obtained the results which are listed in 
Table II, Ill, IV , and IV. The results and source codes using (30) are listed in 
Appendix B for reference purposes. 

Numerical results are obtained for the specific case of L = 400 bits and a = 
4. The levels of diversity considered are M = 1, 2, and 4. A convolutional code of 
rate 1/2, constraint length 7 is used. The values of a, and d are listed in Table I. 


(Ref. 3] 
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Table I. THE COEFFICIENTS 
OF a, AND d 





We will study a ring model for the following cases. 


A. A ONE-RING NETWORK 
A one-ring is used for the network. The user ry locates in a range from R, = 
0 km to R, = 50 km, any other users will vary in the same range with user rp. The 


probability of correct packet capture for one-ring is 


3050 50 


ain) = [f..[A-PQFlrofelrd-Ftnbro-dr,., G9) 
00 O 


where 


Ie 
a O<r, <50 


F,A7,) = 50° 


0 , otherwise 


(40) 


where j = 0, 1, 2, 3, ..., n-1. The calculated values for q(n) are listed in Table IL. 


Table I. SIMULATION RESULTS FOR THE PROBABILITY OF 
CORRECT PACKET CAPTURE q(n) FOR ONE-RING 
NETWORK 





fa fe@ fae fa [as fa fam | ae) 
0587 






M=4 | 1 6318 | 3964 | .2801 | .2132 | .1750 





B. A THREE-RING NETWORK 
The following is an analysis of a three-ring model. Ring one is {0 km < rs 25 
km}; ring two is {25 km < r < 38.5 km} and ring three is {38.5 km < r < 50 km}. 
1. Ring one 


The probability of correct packet capture for ring one is 


where 


20 


74 pe 
USK, < 2) 
feltolRy) = 4 252 ” (42) 


0 , otherwise 


ay iA) eS r,s 50 
flr) = 4 502” (43) 


0 , otherwise 


where j = 1, 2,3 ...n-1. The calculated values of q(n) are listed in Table III. 


Table III. SIMULATION RESULTS FOR THE PROBABILITY OF 
CORRECT PACKET CAPTURE FOR RING ONE 





2. Ring two 


The probability of correct packet capture for ring two is 


50 50 38.5 
gin) = ff f U-Pay ferolR Mar d-Sa nao, (44) 
0 O 25 
where 
ee 1 2385 
FalTolR,) = } 38,52-252° (45) 


0 , otherwise 
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“h Os r, < 50 
fr) =} 57° (46) 


0 ,otherwise 


where j = 1, 2, 3,...n-1. The calculated values of q(n) are listed in Table IV. 


Table IV. SIMULATION RESULTS FOR THE PROBABILITY OF 
CORRECT PACKET CAPTURE FOR RING TWO 


fam fam fae fa [ao [ao [am [ae | 
M=i|i jo fo fo jo [o Jo fo 
.e ee 












3. Ring Three 


The probability of correct packet capture for ring three is 


50 50 50 


g(n) = fu.f [ 1-Pe)fe(rolRaaly)--fxlr-v) (47) 
0 038.5 


v 


where 


2 
a 50 
fdrol) 24 Se ea (48) 


0 , otherwise 


Z2 


a 50 
ears 7, < 
fdr) =4 5? $7 (49) 


0 , otherwise 


where j = 0, 1, 2,...n-1. The calculated values of q(n) are listed in Table V. 


Table V. SIMULATION RESULTS FOR THE PROBABILITY OF 
CORRECT PACKET CAPTURE FOR RING THREE 


——————————— 
aa Tao foo Tea fo Ten Too] 
M=ili [oo [ojo 
M=2 ae On 
m=4 [1 | 2302 | 0210 | 00s Jo fo fo 



















Plots of q(n) versus n for case 1) and each of the rings for case 2) are shown 


in Figures 5, 6, 7, and 8. 
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Figure 5. Probability Of Correct Packet Capture For One-Ring Network 
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Figure 6. Probability Of Correct Packet Capture For Ring One 


25 


San 





_ 
o 
© 
rae) 
A. 
Qt 
2) 
fay 
& 
at aa) 
E 
= 
Pah 
re 
= 
eo) 
E 
ee en ne ee ae 
oO 6 "a (eo joo Teor fe oa 
ainydey yayoeg yOIIIOD JO “Ig 
Figure 7. Probability Of Correct Packet Capture For Ring Two 


26 





0 

Wu 

uM 

Be 

+ 
8) 
~~ 
o 
an 
Sts 
oe 
a) 
ey he 
oO 2 
= 
2 
Z. 
aL 
' S 
t _— 
~ oe 

I 
eS oo Oo Oe SCS 
ginjydey yoxoeg yoot110g jo Aytiqeqoid 
Figure 8. Probability Of Correct Packet Capture For Ring Three 


ws) 


IV. SYSTEM THROUGHPUT FOR A THREE RING ALOHA NETWORK 


The analysis of a ring model for mobile/local radio communications with 
correct packet capture probability has been performed in the previous section. The 
numerical results show that users who are located close to the base station are more 
likely to capture the base station. The throughput for each ring will be calculated by 
using a Markov chain model. The Markov model provides a detailed model of the 
network’s dynamic behavior; therefore, it is a desirable model for the analysis. 

The probability q, that one of i colliding packets will survive a collision and that 
the other i-1 packets will need retransmission later has been calculated in Section IIT. 
The previous system model of an Aloha network with N terminals (users) sharing one 
channel for transmission to one receiver (base station) is used. It is assumed that 
positive acknowledgements from the receiver are never lost, so the user will not 
retransmit unnecessarily. This assumption has been shown to be very realistic since 
modern mobile networks usually have to guarantee a nivty five percent coverage of 
the serving area. The values of q, are calculated in Section III, and it is assumed that 
q, does not equal to zero for i greater than one. This assumption also depends on the 
specific properties of the network. In an undisturbed channel, q, is equal to 1 and qy 


is equal to 0 [Ref. 8]. 
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During any time slot each terminal (user) will be in any of three states: the 
original state O, the transmission state T, or the retransmission state R. The terminal 
will be in state O when there is no packet ready to transmit or if a new packet has 
just arrived while it is waiting for the next time slot. The terminal will be in state T 
when it is actually transmitting a packet, successfully or not. The terminal will be in 
State R after it has transmitted a packet unsuccessfully or when it is waiting for a 
future time slot for retransmission. The terminal will be designated as an O-, T-, or 
R-terminal depending on which state it is in. It 1s assumed that the process which 
generates new packets will be halted in a terminal while it has a packet to transmit. 
The possible state transitions for an individual terminal are shown in Figure 9. 

It is assumed that changes of state can occur at the beginning or at the end of 
the time slot in the Aloha network. An O-terminal can either remain in state O (do 
nothing) or change into a T-terminal at the beginning of a time slot. Upon changing 
into a T-terminal, it starts the transmission of a new packet. The generation rate 
which a terminal changes from state O to state T is denoted py. An R-terminal can 
either stay in state R (do nothing) or change into a T-terminal at the beginning of 
the time slot. In the latter case, it starts the retransmission. The retransmission rate 
when a terminal changes from R to T will be denoted p,. At the end of a time slot, 
a T-terminal can either go to O after a successful transmission or to R after an 


unsuccessful transmission. 


IAS) 


do nothing do nothing 





Figure 9. The Possible State Transition For Each Individual Terminal 


A slotted Aloha network can be described as a homogenous finite Markov 
process with N+1 possible states (0, 1, ...N) corresponding to the number of R- 
terminals at the end of each time slot. 

GivenO<nsNand0<m<N, let 7,,,, denote the state transition probability 
that the network will move from state n to m in one step. There are three different 
types of transitions: 

Type 1, in which the state decreases more than one. This is not possible 
because only one packet will be transmitted successfully in one time slot. 

Type 2, in which the state decreases by one with probability q,. This will happen 
if and only if there is no O-terminal and i R-terminals change into state T for 1 <1 


<n (one packet is successful). 
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Type 3, the state moves to m 2 n. This will happen with probability q,,,.,,, if 
and only if m-n+1 O-terminals and i R-terminals change into state T for 0 <i <n; 
or with probability 1-q,,,,., if and only if m-n O-terminals and 1 R-terminals change 
into state T for 0 <1 < n (no packet successful). 

It is assumed that the terminals are operating independently. The generation 
rate and the retransmission rate are denoted py and p,, respectively. The probability 
of correct packet capture, q,, was obtained in Section HI. The value of 7,,, is 


calculated as [Ref. 8] 








0 ,ifm<n-] 
(pyaar »(" Bpatiep)" |. ¢, , ifm =n-1 
i=1 
oR = Ne oe ca rn a (50) 
i. (|e '. Po D (Fav -P, | 
e ,ifmen 
-m 
: a Py). Giem-n) +> — > Po Gism-n+1 


x 


To calculate the equilibrium state occupation probabilities 7, , 


an arbitrary 


positive number is assigned to 7, . The 7, is calculated as follows [Ref. 6] 








n-] 
Tas = (a. -5 sim (31) 
Tt n-l i=0 | 
Normalizing the 7, , one gets ; 
Th, 
te I (52) 
ri 
i=0 
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The expected throughput of the network in state n (the probability to transmit 


a packet successfully in a time slot) is found to be 


_NA( N-n boy . \Nen-k YO(R i Le, 
f= Ze (vrsng| Por dl 2) DF] 2-4 DY. Gi (53) 


The throughput for a ring is 


N 
$= Son, f, (54) 
n=0 


The system throughput is shown in Figure 10. The curves S1, S2, and S3 are the 
throughputs of ring one, two, and three, respectively. The curve $1+S2+S3 is the 
overall throughput of the network. The source code for Figure 10 is listed in 


Appendix C. 
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V. CONCLUSIONS 


At first, an introduction to Aloha networks have been made. The slotted Aloha 
network is selected in this study, since its system throughput is double the pure Aloha 
system throughput. Then, the probability of correct packet capture of slotted Aloha 
is studied for mobile/local radio communications. An application of slotted Aloha 
in a ring model has been developed. In this model, the collisions do not necessarily 
destroy all the involved packets. With the help of this model, the probability of 
correct packet capture has been investigated. It has been shown that for a ring 
network the system throughput is greater than 1/e, and the ring network is more 
stable under overload. The numerical results in Section HII have been plotted in 
Figures 5, 6, 7, and 8. The probability of correct packet capture of one-ring and 
three-ring networks have been studied. The numerical results and its plot in Figure 
5 show that as the number of users grows from one to six the probability of correct 
packet capture drops rapidly. For each additional user’ , *r six users, the probability 
of correct packet capture does not decrease as quickly as with only a few users. 
When the number of users is large, the system reaches a steady state. The probability 
of correct packet capture increases as M (diversity level) increases. 

Numerical results for a three-ring network show that when the users are located 


in the inner ring (e.g. ring one) they will have a better chance to capture the receiver. 
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As the users move to the outer ring (eg. ring three), the probability that they will 
capture the base station decreases rapidly. The system throughputs for each ring have 
been calculated in Section IV based on the probability of correct packet capture 


obtained previously. 
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APPENDIX A 

As mentioned in Section III, equation (30) is used to calculate q(n) in Tables 
II, HI, IV , and V. The source codes for the results shown in Tables II, III, 1V, and 
V are listed below. Since q(6), q(7), and q(8) require a long computation time, the 
program is written to calculate each value of q(n) separately for large value of n (n 
> 6). 

The source code for the numerical results of Table II is listed below. The 
number of users (nuser) is initialized to 2, 3,..., and 8 to calculate q(2), q(3), ..., q(8), 
respectively. In the listed program the maximum number of users (nuser) is equal to 


eight. For a one ring network, R, and R, are varied from 0 km to 50 km. 
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%(* this program calculates q(8) for one ring network *) 
%(* the coefficient values for d and a, *) 

d =(10 12 14 16 18 10}; 

ad={11 38 193 1331 7275 40406); 


%(* the number of interval for integration *) 
interv=5; 


%(* the values for R, and R, *) 
r1=[0 0}; 
r2=[50 50}; 


%Yo(* set value of a, L, number of user, and m1 for M *) 
alpha=4; 

L=400; 

nuser = 8; 

m1=[{1 2 4]; 


%(* calculate the grid points of integration *) 
delta(1)=(r2(1)-r1(1))/interv; 
delta(2) =(1r2(2)-r1(2))/interv; 
r(1,1)=1r1(1)+delta(1)/2; 
r(2,1) =r1(2) + delta(2)/2; 
for i=2:interv 

r(1,1) =r(1,i-1) + delta(1); 
end | 
for 1=2:interv 

r(2,i) =1r(2,i-1) + delta(2); 
end 


%(* calculate the probability density function fp equation (7) *) 
for i= l:interv | 
fr(1,1) =2* r(1,1)/(r2(1)*2-r1(1)*2); 
en 
for 1= 1:interv 
fr(2,1) =2*r(2,1)/(r2(1)*2-r1(2)*2); 
end 
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%Yo(* this loop calculate q(8) for all 
Yovalue of M=1, 2, and 4, equation (38) *) 
for i12=1:3 : 
m=m1(i2); 
sum = 0; 
for a=l:interv 
for b= 1:interv 
for c=1:interv 
for f=1:interv 
for g=1:interv 
for h=1:interv 
for k=1:interv 
for o=l:interv 
index=[a bc fg hk o]; 
muy =0; 


%(* calculate value for u, equation (27) *) 
for 11 =1:nuser-1 
muy = muy + (r(1,index(1)/r(2,index(i1 + 1)))*alpha; 
end 
muy = 1/sqrt(1+ muy); 


%(* calculate value for P,, equation (25) *) 
pb =(1-muy)/2 


%o(* calculate the value for P,.,, equation (26) *) 
pbm=0; 
for cnt=0:m-1 
pbm=pbm + (comb((m-1+ cnt),cnt)*(((1 + muy)/2)*cnt)); 
end 
pbm=(pb*m)*pbm; 


%(* calculate the value for P, for d even, equqation (36) *) 
pd= (2*sqrt(pbm*(1-pbm))).“d; 


%(* calculate P., equation (29) *) 


pe=ad*pd’; 
if pe>1,pe=1;,end 
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%(* calculate the integration for q(n), equation (38) *) 
inc=(1-pe)*L* fr(1,0)*fr(2,a) *fr(2,b) *fr(2,c) *fr(2,f); 
inc=inc* fr(2,g)*fr(2,h)*fr(2,k); 
sum =sum + inc; 
end 
end 
end 
end 
end 
end 
end 
end 
sum =sum*delta(2)*(nuser-1)*delta(1); 
end 


F 

The source code for the numerical results of Table III is listed below. For ring 
one, R, varies from 0 km to 25 km, and R, varies from 0 km to 50 km. This source 
code is used to calculate q(8). To compute q(2), q(3),....q(n), the number of users 


(nuser) is equal to 2, 3,..., n, respectively. 


oy 


%(* this program calculates q(8) for ring one *) 
%o(* the coefficient values for d and a, *) 

d =[10 12 14 16 18 10]; i 
ad=[11 38 193 1331 7275 40406]; 


%(* the number of interval for integration *) 
interv=5; 


Yo(* the values for R, and R, *) 
r1=[0 0]; 
r2=[25 50]; 


Jo(* set value of a, L, number of user, and m1 for M *) 
alpha =4; 

L=400; 

nuser=8; 

mi=[1 2 4]; 


Yo(* calculate the grid points of integration *) 
delta(1) =(r2(1)-r1(1))/interv; 
delta(2) =(r2(2)-r1(2)) /interv; 
r(1,1)=r1(1)+delta(1)/2; 
(2,1) =r1(2) + delta(2)/2; 
for i=2:interv 

r(1,i) =r(1,i-1) + delta(1); 
end 
for i=2:interv 

r(2,1) =1r(2,i-1) + delta(2); 
end 


%(* calculate the probability density function fg eq. ‘.ion (7) *) 
for 1=1:interv 
fr(1,i) =2*r(1,i)/(r2(1)*2-r1(1)*2); 
end 
for i= 1:interv 
fr(2,1) =2*r(2,1)/(r2(1)°2-r1(2)*2); 
end 
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%(* this loop calculates q(8) for all 
JYovalue of M=1, 2, and 4, equation (38) *) 
fon i2= 1:3 
m=m1(12); 
sum =0; 
for a= 1:interv 
for b=1:interv 
for c= l:interv 
for f= l:interv 
for g=l:interv 
for h=1:interv 
for k=1:1nterv 
for o= ltinterv 
index=[a bcfghk o]; 
muy = 0; 


%(* calculate value for 4, equation (27) *) 
for 11=1:nuser-1 
muy = muy + (r(1,index(1)/r(2,index(i1 + 1)))*alpha; 
end 
muy = 1/sqrt(1+ muy); 


%o(* calculate value for P,, equation (25) *) 
pb=(1-muy)/2 


%(* calculate the value for P,,,, equation (26) *) 
pbm=0; 
for cnt=0:m-1 
pbm=pbm-+(comb((m-1+cnt),cnt)*(((1+ muy)/2)“ent)); 
end 
pbm =(pb*m)*pbm; 


%(* calculate the value for P, for d even, equqation (36) *) 
pd =(2*sqrt(pbm*(1-pbm))).“d; 


%(* calculate P,, equation (29) *) 


pe=ad*pd’; 
if pe>1,pe=1;,end 
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%(* calculate the integration for q(n), equation (38) *) 
inc=(1-pe)*L*fr(1,0)*fr(2,a) *fr(2,b)* fr(2,c)* fr(2,f); 
inc= inc*fr(2,g)*fr(2,h)*fr(2,k); 
sum =sum+ inc; 
end 
end 
end 
end 
end 
end 
end 
end 
sum =sum*delta(2)*(nuser-1)*delta(1); 
end 


LESSEE —EEL— EE eS EEE Sg SSE aa 

The source code for the numerical results of Table IV is listed below. For ring 
two, R, varies from 25 km to 38.5 km, and R, varies from 0 km to 50 km. This 
program is used to calculate q(8). To calculate q(2), q(3),.... q(n), the number of 


users (nuser) is equal to 2, 3,..., n, respectively. 
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%o(* this program will calculate q(8) for ring two *) 
%o(* the coefficient values for d and a, *) 

d =[10 12 14 16 18 10]; 

ad=[11 38 193 1331 7275 40406]; 


Yo(* the number of interval for integration *) 
interv=5; 


Yo(* the values for R, and R, *) 
r1=[25 0]; 
r2 =(38.5 50]; 


Yo(* set value of a, L, number of user, and m1 for M *) 
alpha=4; 

L=400; 

nuser=8; 

ml=[1 2 4]; 


%o(* calculate the grid points of integration *) 
delta(1)=(r2(1)-r1(1))/interv; 
delta(2) =(1r2(2)-r1(2))/interv; 
r(1,1)=r1(1)+delta(1)/2; 
r(2,1)=r1(2)+delta(2)/2; 
for i=2:interv 

r(1,i)=r(1,i-1) + delta(1); 
end 
for 1=2:interv 

r(2,1) =1(2,i-1) + delta(2); 
end 


%o(* calculate the probability density function fp equation (7) *) 
for i=1:interv 
fr(1,i) =2*r(1,i)/(r2(1)*2-r1(1)*2); 
end 
for 1=1:interv 
fr(2,1) =2*r(2,1)/(12(1)*2-r1(2)*2); 
end 
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%(* this loop calculates q(8) for all 
Yovalue of M=1, 2, and 4, equation (38) *) 
for 12=1:3 
m=m1(i2); 
sum =Q; 
for a=1:interv 
for b=1:interv 
for c= l:interv 
for f=1:interv 
for g=1:interv 
for h=1:interv 
for k= 1:interv 
for o= 1:interv 
index=[a bc f ghk o]; 
muy =0; 


%(* calculate value for y, equation (27) *) 
for 11 =1:nuser-1 
muy = muy + (r(1,index(1)/r(2,index(i1 + 1)))*alpha; 
end 
muy = 1/sqrt(1+ muy); 


Yo(* calculate value for P,, equation (25) *) 
pb =(1-muy)/2 


Yo(* calculate the value for P,,,, equation (26) *) 
pbm=0; 
for cnt=0:m-1 
pbm= him+ (comb(a Tee ent)*(((1+ muy)/2)*cnt)); 
end 
pbm=(pb"m)*pbm; 


%(* calculate the value for P, for d even, equqation (36) *) 
pd =(2*sqrt(pbm*(1-pbm))).“d; 


Yo(* calculate P,, equation (29) *) 


pe=ad*pd’; 
if pe > 1,pe=1;,end 
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%(* calculate the integration for q(n), equation (38) *) 
inc=(1-pe)*L* fr(1,0)* fr(2,a) *fr(2,b) * fr(2,c)* fr(2,f); 
inc=inc*fr(2,g)* fr(2,h) *fr(2,k); 
sum= sum-+inc; 
end 
end 
end 
end 
end 
end 
end 
end 
sum =sum*delta(2)*(nuser-1)*delta(1); 
end 





The source code for the numerical results of Table V is listed below. For ring 
three, R, varies from 38.5 km to 50 km, and R, varies from 0 km to 50 km. This 
program is used to calculate q(8). To calculate q(2), q(3),.... q(n), the number of 


users (nuser) is equal to 2, 3,..., n, respectively. 
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%(* this program will calculate q(8) for ring three *) 
Yo(* the coefficient values for d and a, *) 

d =[10 12 14 16 18 10]; 

ad=[11 38 193 1331 7275 40406]; 


Yo(* the number of interval for integration *) 
interv=5; 


Yo(* the values for R, and R, *) 
r1=[38.5 0]; 
r2=(50 SO]; 


Yo(* set value of a, L, number of user, and m1 for M *) 
alpha =4; 

L=400; 

nuser = 8; 

m1=[1 2 4]; 


Yo(* calculate the grid points of integration *) 
delta(1) =(r2(1)-r1(1))/interv; 
delta(2) = (r2(2)-r1(2)) /interv; 
r(1,1)=r1(1) + delta(1)/2; 
r(2,1)=1r1(2) + delta(2)/2; 
for i=2:interv 

r(1,1) =r(1,1-1) + delta(1); 
end 
for i=2:interv 

r(2,1) =r(2,i-1) + delta(2); 
end 


Y%o(* calculate the probability density function fp eq ” ion (7) *) 
for 1= l:interv 
fr(1,i) =2*r(1,1)/(r2(1)*2-r1(1)*2); 
end 
for 1=1:interv 
fr(2,1) =2*1r(2,1)/(12(1)°2-r1(2)*2); 
end 
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%(* this loop calculates q(8) for all 
Yovalue of M=1, 2, and 4, equation (38) *) 
rOimaz — 1:3 
m=m1(i2); 
Sum = 0; 
for a=l1:interv 
for b= l:nterv 
for c= l:interv 
for f=1:interv 
for g=l:interv 
for h=1:mterv 
for k= 1:interv 
for o= l:interv 
index=[a bc fghk o]; 
muy =0; 


%(* calculate value for u, equation (27) *) 
for 11 = 1:nuser-1 
muy = muy + (r(1,index(1)/r(2,index(i1 + 1)))*alpha; 
end 
muy = 1/sqrt(1+muy); 


%o(* calculate value for P,, equation (25) *) 
pb =(1-muy)/2 


%o(* calculate the value for P,,,, equation (26) *) 
pbm =0; 
for cnt=0:m-1 
pbm =pbm+ (comb((m-1 + cnt),cnt)*(((1 + muy)/2)*ent)); 
end 
pbm =(pb*m)*pbm; 


%(* calculate the value for P, for d even, equqation (36) *) 
pd =(2*sqrt(pbm*(1-pbm)))"d; 


%(* calculate P,, equation (29) *) 


pe=ad*pd’; 
if pe >1,pe=1;,end 
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a SS EE SS E—EEEE—E——EEEEET 
7o(* calculate the integration for q(n), equation (38) *) 
inc=(1-pe)*L* fr(1,0)*fr(2,a)*fr(2, Dae c)*fr(2,f); 
inc= inc*fr(2,g)* fr(2,h)* fr(2,k); | 
sum =sum + inc; 

end 
end 
end 
end 
end 
end 
end 
end 
sum =sum*delta(2)*(nuser-1)*delta(1); 
end 
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APPENDIX B 


Equation (36) is used to calculate q(n) since computing q(n) for n 2 6 from 
equation (30) is prohibitively time-consuming. The following programs use (30) for 
P, to calculate q(n). The following programs are listed for reference purposes. The 
source codes and the results are listed below. Since q(n) for n > 6 takes a long time 
to calculate, the program is written to calculate only up to q(6). 

The number of users (nuser) is initialized to 2, 3,..., 6 to calculate q(2), q(3), 
..., q(6), respectively. In the listed program the number of users (nuser) is equal to 
six. For a one-ring network, R, and R, vary from 0km to 50km. The source code and 


its simulation results are listed below and in Table VI, respectively. 
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Y(* this program calculates q(6) for one ring network *) 
%(* the coefficient values for d and a, *) 

d =[10 12 14 16 18 10}; 

ad=[11 38 193 1331 7275 40406); 


Yo(* the number of interval for integration *) 
interv=5; 


Yo(* the values for R, and R, *) 
r1=[(0 0); 
r2=[50 50]; 


Yo(* set value of a, L, number of user, and m1 for M *) 
alpha =4; 

L=400; 

nuser =8; 

ml=[1 2 4]; 


%(* calculate the grid points of integration *) 
delta(1) =(r2(1)-r1(1))/interv; 
delta(2) =(r2(2)-r1(2))/interv; 
r(1,1)=r1(1)+delta(1)/2; 
r(2,1) =r1(2)+delta(2)/2; 
for i=2:interv 

r(1,i1) =r(1,i-1) + delta(1); 
end 
for 1=2:interv 

r(2,1) =r(2,i-1) + delta(2); 
end 


Yo(* calculate the probability density function fp equation (7) *) 
for i=1:interv 
fr(1,1) =2*r(1,i)/(r2(1)*2-r1(1)*2); 
end 
for i=1:interv 
fr(2,1) =2*r(2,1)/(12(1)*2-r1(2)*2); 
end 
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%(* this loop calculate q(6) for all 
Yovalue of M=1, 2, and 4, equation (38) *) 
hor 12.= 1:3 
m=m1(i2); 
sum =0; 
for a=l:interv 
for b= 1:interv 
for c=1:interv 
for f=1:nterv 
for g=1:interv 
for o= l:interv 
index=[a b c f g oJ; 
muy =0; 


%(* calculate value for yu, equation (27) *) 
for 11=1:nuser-1 
muy = muy + (r(1,index(1)/r(2,index(i1 + 1)))*alpha; 
end 
muy = 1/sqrt(1+ muy); 


%(* calculate value for P,, equation (25) *) 
pb=(1-muy)/2 


%(* calculate the value for P,,,, equation (26) *) 
pbm =0; 
for cnt=0:m-1 
pbm =pbm+(comb((m-1+cnt),cnt)*(((1+ muy)/2)*ent)); 
end 
pbm =(pb*m)*pbm; 


%(* calculate the value for P; for d even, equqation (30) *) 
for q=1:6 
pd(q) =0.5*comb(d(q),d(q)/2)*pbm*(d(q)/2); 
pd(q) =pd(q)*(1-pbm)*(d(q)/2); 
for e=1+d(q)/2:d(q) 
Pala) = Paka) + (comb(a¢q),e)"{(pbm ¢)"C-pom) (alae) 
en 
end 


5) II 


%(* calculate P., equation (29) *) 
pe=ad*pd’; 
if pe>1, pe=1;, end 


%o(* calculate the integration for q(n), equation (38) *) 
inc=(1-pe)*L* fr(1,0)*fr(2,a)*fr(2,b) *fr(2,c) *fr(2,f); 
inc=inc*fr(2,g); 

sum =sum+ inc; 

end 
end 
end 
end 
end 
end 
sum =sum*delta(2)*(nuser-1)*delta(1); 
end 





Table VI. SIMULATION RESULTS FOR THE PROBABILITY OF 
CORRECT PACKET CAPTURE FOR ONE-RING NETWORK 





The source code and the simulation results are listed below and in Table VU, 
respectively. For ring one R, is varied from 0 km to 25 km, and R, is varied from 0 
km to 50 km. This source code is used to calculate q(6). To compute q(2), 


q(3),....q(n), the number of users (nuser) is 2,3,....n respectively. 


SZ 





Y%o(* this program calculates q(6) for ring one *) 
Yo(* the coefficient values for d and a, *) 

d =[10 12 14 16 18 10]; 

ad=[11 38 193 1331 7275 40406]; 


%o(* the number of interval for integration *) 
interv=S; 


Yo(* the values for R, and R, *) 
r1=[0 0]; 
r2=[(25 SO}; 


Yo(* set value of a, L, number of user, and m1 for M *) 
alpha=4; 

L=400; 

nuser =8; 

m1=[1 2 4]; 


%(* calculate the grid points of integration *) 
delta(1) = (r2(1)-r1(1)) /interv; 
delta(2) =(r2(2)-r1(2))/interv; 
r(1,1)=1r1(1)+ delta(1)/2; 
r(2,1)=r1(2) + delta(2)/2; 
for i=2:interv 

r(1,i)=r(1,i-1) + delta(1); 
end 
for 1=2:interv 

r(2,i)=1(2,i-1) + delta(2); 
en 


%(* calculate the probability density function fp equation (7) *) 
for i=1:interv 
fr(1,i) =2*r(1,i)/(r2(1)*2-r1(1)*2); 
end 
for i=1:interv 
fr(2,1) =2*1(2,i) /(r2(1)*2-r1(2)*2); 
en 
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%(* this loop calculates q(6) for all 
Yovalue of M=1, 2, and 4, equation (38) *) 
for 12=1:3 
m=m1(i2); 
sum=0; 
for a= l:interv 
for b= 1:interv 
for c=1:interv 
for f=1:interv 
for g=1:interv 
for o=1:interv 
index=[a b cf g oJ; 
muy =0; 


%(* calculate value for p, equation (27) *) 
for i1=1:nuser-1 
muy = muy + (r(1,index(1)/r(2,index(il + 1)))“alpha; 
end 
muy = 1/sqrt(1+ muy); 


%(* calculate value for P,, equation (25) *) 
pb =(1-muy)/2 


%(* calculate the value for P,,,, equation (26) *) 
pbm=0; 
for cnt=0:m-1 
pbm =pbm + (comb((m-1+ cnt),cnt)*(((1 + muy) /2)“cnt)); 
end 
pbm =(pb*m)*pbm; 


%(* calculate the value for P, for d even, equqation (30) *) 
for q=1:6 
pd(q) = 0.5*comb(d(q),d(q)/2)*pbm™(d(q)/2); 
pd(q)=pd(q)*(1-pbm)*(d(q)/2), 
for e=1+d(q)/2:d(q) 
Pata) Pal) + (comb ee 1-pbm)*(d(q)-e)); 
en 
end 
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%(* calculate P., equation (29) *) 
pe=ad*pd’; 
if pe>1, pe=1;, end 


%(* calculate the integration for q(n), equation (38) *) 
inc=(1-pe)*L* fr(1,0)*fr(2,a)*fr(2,b)* fr(2,c) *fr(2,f); 
inc=inc*fr(2,g); 

sum =sum + inc; 

end 
end 
end 
end 
end 
end 
sum =sum*delta(2)*(nuser-1)* delta(1); 
end 





Table VII. SIMULATION RESULTS FOR THE PROBABILITY OF 
CORRECT PACKET CAPTURE FOR RING ONE 





The source code and the simulation results are listed below and in Table VIII, 
respectively. For ring two, R, is varied from 25 km to 38.5 km, and R, is varied from 
0 km to 50 km. This program is used to calculate q(6). To calculate q(2), q(3),..., 


q(n), the number of users (nuser) is 2, 3,..., n respectively. 


> 


%(* this program will calculate q(6) for ring two *) 
Ge the coefficient values for d and a, 2. 

d =[10 12 14 16 18 10}; 
ad=[11 38 193 1331 7275 40406]; 


Yo(* the number of interval for integration *) 
interv=5; 


Y(* the values for R, and R, *) 
=[25 0]; 
r2=[38.5 50]; 


%(* set value of a, L, number of user, and m1 for M *) 
alpha=4; 

L=400; 

nuser =8; 

ml=[1 2 4]; 


%o(* calculate the grid points of integration *) 
delta(1)=(r2(1)-r1(1))/interv; 
delta(2) =(12(2)-r1(2))/intery; 
r(1,1)=1r1(1)+delta(1)/2; 
r(2,1) =r1(2) + delta(2)/2; 
for i=2:interv 

r(1,i) =r(1,i-1) + delta(1); 
end 
for 1=2:interv 

r(2,1) =r(2,i-1) + delta(2); 
end 


%(* calculate the probability density function fp equation (7) *) 
for i=1:interv 
fr(1,1) =2*r(1,i)/(r2(1)*2-11(1)*2); 
en 
for 1= 1:interv 
fr(2,1) =2*1r(2,i) /(r2(1)*2-r1(2)*2); 
end 
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%(* this loop calculates q(6) for all 
Yovalue of M=1, 2, and 4, equation (38) *) 
for i2=1:3 
m=m1(i2); 
sum=0; 
for a= l:interv 
for b= 1:interv 
for c=1:interv 
for f=1:interv 
for g=1:interv 
for o=1:interv 
index=[a b cf g oJ; 
muy =0; 


%(* calculate value for u, equation (27) *) 
for i11=1:nuser-1 
muy = muy + (r(1,index(1)/r(2,index(i1 + 1)))*alpha; 
end 
muy = 1/sqrt(1+ muy); 


%(* calculate value for P,, equation (25) *) 
pb=(1-muy)/2 


%(* calculate the value for P,,,, equation (26) *) 
pbm=0; 
for cnt=0:m-1 
pbm =pbm + (comb((m-1+ ent),cnt)*(((1 + muy)/2)*cnt)); 
end 
pbm =(pb"m)*pbm; 


%(* calculate the value for P; for d even, equqation (30) *) 
for q=1:6 
pd(q) =0.5*comb(d(q),d(q)/2)*pbm(d(q)/2); 
pd(q) =pd(q)*(1-pbm)*(d(q)/2); 
for e=1+d(q)/2:d(q) 
BORD serena + Koomibccke).c) (poem 2 Pom) Aaa): 
en 
end 


Dy. 


%(* calculate P., equation es) A 
pe=ad*pd’; 
if pe>1, pe=1;, end 


%(* calculate the integration for q(n), equation (38) *) 
inc =(1-pe)*L* fr(1,0)*fr(2,a)*fr(2,b)*fr(2,c) *fr(2,f); - 
inc=inc*fr(2,g); | 

sum =sun + inc; 





end 
end 
end 
end 
end 
end 
sum =sum*delta(2)*(nuser-1)*delta(1); 
end 
Table VIII. SIMULATION RESULTS FOR THE PROBABILITY OF 


CORRECT PACKET CAPTURE FOR RING TWO 


Pee fact) BCom =a 
ee ee 
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The source code and the simulation results are listed below and in Table IX, 
respectively. For ring three, R, is varied from 38.5 km to 50 km, and R, is varied 
from 0 km to 50 km. This program is used to calculate q(6). To calculate q(2), 


q(3),.... q(n), the number of users (nuser) is 2, 3,..., mn respectively. 
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%(* this program will calculate q(6) for ring three *) 
Yo(* the coefficient values for d and a, *) 

d =(10 12 14 16 18 10]; 

ad=[11 38 193 1331 7275 40406); 


%(* the number of interval for integration *) 
interv=5; 


%(* the values for R, and R, *) 
r1=(38.5 0]; 
r2=(50 50]; 


Yo(* set value of a, L, number of user, and m1 for M *) 
alpha=4; 

L=400; 

nuser =8; 

m1=[1 2 4]; 


Yo(* calculate the grid points of integration *) 
delta(1) =(r2(1)-r1(1)) /intery; 
delta(2) =(r2(2)-r1(2))/interv; 
r(1,1)=r1(1)+delta(1)/2; 
r(2,1)=11(2) + delta(2)/2; 
for i=2:interv 

r(1,1) =1(1,i-1) + delta(1); 
end 
for i=2:interv 

1(2,i) =1(2,i-1) + delta(2); 
end 


. Y%(* calculate the probability density function fp equation (7) *) 
for i= L:interv 
fr(1,i1) =2*r(1,i)/(r2(1)*2-r1(1)*2); 
end 
for i=1:interv 
fr(2,1) =2*r(2,1)/(12(1)*2-r1(2)*2); 
end 
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%(* this loop calculates q(6) for all 
Yovalue of M=1, 2, and 4, equation (38) *) 
for 12=1:3 
m=m1(i2); 
sum=0Q; 
for a=l:interv 
for b=1:interv 
for c= l:interv 
for f=1:interv 
for g=1:interv 
for o= l:interv 
index=[a b c f g o]; 
muy =0; 


%(* calculate value for yu, equation (27) *) 
for 11=1:nuser-1 
muy =muy + (r(1,index(1)/r(2,index(il + 1)))“alpha; 
end 
muy = 1/sqrt(1+ muy); 


%(* calculate value for P,, equation (25) *) 
pb=(1-muy)/2 


%o(* calculate the value for P,.,, equation (26) *) 
pbm=0; 
for cnt=0:m-1 
pbm=pbm + (comb((m-1+ ent),cnt)*(((1 + muy) /2)*ent)); 
end 
pbm=(pb*m)*pbm; 


%(* calculate the value for P, for d even, equaation (30) *) 
for q=1:6 
pd(q) =0.5*comb(d(q),d(q)/2)*pbm”(d(q)/2); 
pd(q) =pd(q)*(1-pbm)*(d(q)/2); 
for e=1+d(q)/2:d(q) 
pay = pd(q) + (comb(d(q),e)*(pbm*e)*(1-pbm)"(d(q)-e)); 
en 
end 
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%(* calculate P,, equation (29) *) 
pe=ad*pd’; 
if pe>1, pe=1;, end 


%(* calculate the integration for q(n), equation (38) *) 
inc=(1-pe)*L*fr(1,0)*fr(2,a) *fr(2,b) * fr(2,c)* fr(2,f); 
inc=inc*fr(2,g); 

sum = sum + inc; 

end 
end 
end 
end 
end 
end 
sum =sum* delta(2)*(nuser-1)*delta(1); 
end 





Table IX. SIMULATION RESULTS FOR THE PROBABILITY OF 
CORRECT PACKET CAPTURE FOR RING THREE 





2518 0447 





Plots of q(n) versus n for one-ring and three-ring systems are shown in Figure 


ee 12, 13, and 14. 
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Figure 11. Probability Of Correct Packet Capture For One-Ring Network 
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Figure 12. Probability Of Correct Packet Capture For Ring One 
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Figure 13. Probability Of Correct Packet Capture For Ring Two 
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Figure 14. Probability Of Correct Packet Capture For Ring Three 
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APPENDIX C 


SS SS SSeS SS LE OEE LL SS ——E———— 
Yo(* initialization *) 
nuser = 15; 
p0=0.1; 
pr =0.002; 
q=[0 1 0.9042 0.8099 0.7197 0.6350 0.5580 0.4901 0.4615 0.4285]; 


Yo(* steady state value of q(n) *) 
for 1=11:nuser+4 
q(i) = 0.3852; 
end 
%o(* set N=nuser, calculate 7, ,,, it is called pi(n,m) *) 
N =nuser; 
for n=0:N 
nl=n+1; 
for m=0:N 
mi=m+1; 
if m<n-1,pi(ni,m1) =0; 


n,m? 


elseif m= =n-1 
sum1 =0; 
for i=0:n 
11=1+1; 
suml =sum1+comb(n,i)*pr“i*(1-pr)*(n-1)*qQ@i1); 
end 
pi(nl,m1)=(1-p0)*(N-n) *sum1; 


elseif m>=n 

sum1 =0; 

temp1=comb(N-n,N-m)*(1-p0)*(N-m-1)*p0°(m-n); 

for i=0:n 
il=i+1 
term =comb(n,1)*(1-pr)*(n-i)* pr; 
term=term*((1-p0)*(1-q(i1+ m-n)) + (N-m)*p0*q(i1+m-n+ 1)/(m-n+ 1)); 
sum1l=sum1+term; 

end 

pi(nl,m1)=temp1 *sum1; 

SSS —  —eeeeeeeeeeeee———————EEE EE  — —————————— ee 


66 


%(* calculate 7, , it is called pis(n) *) 
for n=1:N 
ni=n+1; 
Sel. — Ue 
for 1=0:n-1 
1l=i+1; 
sum1=sum1 + (pis(i1)*pi(ii,n1-1)); 
end 
temp =pis(n1-1)-sum1; 
pis(n1) =(1/pi(n1,n1-1))*temp; 
if pis(n1) <0,pis(n1) =0;,end 
end 


%(* calculate 7, it is called pin *) 
den=sum(pis); 
pin=pis/den; 
%o(* calculate f, *) 
for n=0:N 
nl=n+1; 
temp =0; 
total =0; 
for k=0:N-n 
k1=k+1; 
sum1 =0; 
for 1=0:n 
sum1=sum1+(comb(n,i)*pr*i*(1-pr)“(n-i)*q(it+k1); 
end 
temp =comb(N-n,N-n-k)*p0°k(1-p0)°(N-n-k); 
term=sum1*temp; 
total =total+ term; 
end 
f(n1) =total; 
end 
%(* calculate s *) 
sum1=0; 
for 1=0:N 
il=i+1; 
sum1=sum1 + (pin(i1)*f(i1)); 
s(i1) =sum1; 


end — SSSSS=S=S—0@MM0 SSeS 
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%o(* calculate the Comb(n,m) *) 
function y=comb(n,m) 
y=xfact(n+ 1)/(xfact(m + 1)*xfact(n-m +!)); 


%o(* calculate the factorial, it is call fact(n) *) 
function y=fact(x) 
ifx== 

y= 
else 

Wek 

for i=1:x 

y=y"l, 

end 

end 


Y(* the variable xfact *) 
x=0:52 
for i=1:53 

xfact(i) =fact(x(i)); 
end 


global xfact 
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